function  answer = ejer3()

	x = 1:1:5;
	f = 1./x;
	y = [];
	
	for i = 1 : length(x)
		y = [y, f(x(i))];
	endfor
	
	
	z = [];
	for j = 1:101
		z = [z,24/25 + j/25];
	endfor
	
	w = newton(x,y,z);
	
	plot(z,w,'g');
	hold on 
	plot(x,f);
  
endfunction 


